Non-Uniform Cosmological Birefringence and Active Galactic Nuclei 
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Cosmological birefringence, a rotation by an angle a of the polarization of photons as they prop- 
agate over cosmological distances, is constrained by the cosmic microwave background (CMB) to 
be \a\ < 1° (la) out to redshifts z ~ 1100 for a rotation that is uniform across the sky. How- 
ever, the rotation angle a(9, <f>) may vary as a function of position (0, <j>) on the sky. Here I discuss 
how a position-dependent rotation can be sought in current and future AGN data. An upper limit 

^a 2 ^ 1//2 < 3.7° to the scatter in the position-angle-polarization offsets in a sample of only N = 9 



AGN already constrains the rotation spherical-harmonic coefficients to (47r) 



-1/2 



< 3.7° and 



constrains the power spectrum for a in modeis where it is a stochastic field. Future constraints can 
be improved with more sources and by analyzing well-mapped sources with a tensor-harmonic de- 
composition of the polarization analogous to that used in CMB pofarization and weak gravitational 
lensing. 



PACS numbers: 98.80.-k 

Introduction. There is a very active quest to un- 
derstand dark energy |l|, and quintessence models 
[H provide a promising set of effective theories. A 
pseudo-Nambu-Goldstone field provides an attractive 
quintessence candidate, and such a field should have a 
coupling to the Chern-Simons term of electromagnetism 
0| . This coupling gives rise to cosmological birefringence 
(CB), a frequency-independent rotation by some angle 
a, of the linear polarization of photons as they propa- 
gate over cosmological distances 0, 0] ■ I n the simplest 
models, the rotation angle a is uniform across the sky, 
in which case CB gives rise to parity-violating TB and 
EB correlations in cosmic microwave background (CMB) 
maps Null searches for such correlations now con- 
strain the rotation over the redshift range < z < 1100 
to be a = -0.25° ±0.58° 0. 

However, several recent papers have introduced 
quintessence models in which the rotation angle a( 
varies as a function of position (9, <j>) on the sky 
and similar phenomena may arise in some dark-matter 
models Q. Refs. @ have described how to measure 
this rotation angle, as a function of position on the sky, 
with the CMB, but the algorithm has not yet been ap- 
plied to data. The WMAP satellite should be sensi- 
tive (at la) to spherical-harmonic coefficients of the ro- 



tation as small as (Alt) 



-1/2 



Oilr, 



2.3°, for I < 400, 



and the recently launched Planck satellite should reach 
(4 7 r)- 1 / 2 |Q:i m | - 0.07° for I < 800 |]. 

Radio 

SEMI and UV 0,13 data on active galactic 
nuclei (AGN) can also be used to search for CB. AGN are 
often elongated and polarized. While AGN may be com- 
plicated objects, symmetry considerations suggest that 
on average the mean polarization should be aligned or 
perpendicular to the position angle of the source. CB 
would, by rotating the polarization, give rise to a nonzero 
mean offset between the position angles and polarizations 
measured in a large number of sources 0, E3, EH ■ Like- 



wise, if a detailed map of the intensity and polarization 
of an individual source can be made 12J, Il3l . Il5j , then on 
average, the intensity gradients and polarization within 
that source should be aligned or perpendicular, if there 
is no CB. This more detailed analysis may allow compet- 
itive, or even stronger, constraints to a from a smaller 
number of sources. 

One precisely imaged radio source (3C 9) at a redshift 
z ~ 2 constrains a — 2° ± 3° out to this distance IH 17|. 
A stronger bound, a = —0.6° ± 1.5°, can be obtained 
from a larger number of well-mapped radio sources, but 
only at smaller redshifts [3[. A recent UV sample [16( of 
eight AGN at redshifts z > 2 constrains a = -0.7°±2.0°. 
In this Letter, I show that AGN can be used to constrain 
the multipole moments ai m for a non-uniform rotation. 

There was a brief flurry in the 1990s of searches for 
a rotation with a dipole dependence on position [ill — 
14 1, following a claimed detection [lfj]. Here I revisit 



and update such measurements and generalize to higher- 
l moments. I search the recent UV da ta Jig , combined 
with the radio constraint from 3C 9 [l3j . and find no 
evidence for any rotation with a dipole or quadrupole 
dependence on position. I constrain the ai m (for any 
I) to (47r) _1 / 2 a; m < 3.7° (ler), and I place a constraint 
to the power spectrum for a for theories that predict 
that a(9, </>) is a stochastic field. To preface, I discuss 
the derivation of the constraint, a — —0.7° ± 2.0° (to 
redshifts z ~ 2), from the UV data, a result that is 
strengthened to a = —0.1° ± 1.7° if the radio data on 
3C 9 is included. Finally, I discuss how the analysis of 
high-resolution intensity-polarization maps of individual 
sources can be optimized, using tensor-harmonic tech- 
niques similar to those for CMB polarization and weak 
gravitational lensing. 

Prelude: A uniform rotation. Table |T] reproduces data 
on 8 UV sources from Ref. [lij] as well as radio data 
on 3C 9 [H \vf\ . Listed there are the positions (6i,4>i), 
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TABLE I: The 6i-(f>t coordinates, offsets Ui and measurement 
errors ai for the eight sources listed in Ref. pi| plus (the last) 
the radio source 3C 9 (from Ref. [l3|). 



position-angle-polarization offsets a.;, and measurement 
errors Oi to these offsets. 

Let us first test with this data whether there is a rota- 
tion, by an angle a, that is uniform across the sky. We 
will also determine the scatter o~ p in the measurements 
of a. To estimate the mean offset from the data, we use 
the minimum- variance estimator, 



(i) 



The error to our measurement of a is then inverse root 
of the denominator. We find for the UV data in Table U 
a = — 0.7 ±2.0. Including the radio source 3C 9 improves 
the minimum- variance constraint to a — —0.1° ± 1.7°. 
In general, it could be that the sample contains a sig- 
nificant intrinsic scatter in the offsets, in which case the 
minimum-variance error would underestimate the true 
error in a. For example, one extremely well-measured 
and statistically-significant nonzero offset might suggest 
nonzero CB, but could alternatively be due to an intrin- 
sic offset in the source. Still, the measured offsets for the 
current sample are well within their measurement errors, 
and so the error obtained here is probably sound. The 
error cr 2 = J2i( a i ~ oi) 2 /[N(N — 1)] obtained from the 
measured dispersion is in fact a bit smaller, o~ a = 1.4°, 
suggesting that the reported measurement errors in this 
sample may be a bit high and the true constraint a bit 
stronger. 

The scatter a v in the offsets is determined from the 
data via, 



E 



-r 4 /E-r 2 



(2) 

The 9 sources in Table U result in a p — 2.9°, a result 
that will be used below. Note that the weighted estimate 
of the scatter in Eq. ([2]) is a bit smaller than the value 
4.4° obtained if an unweighted estimator, <r 2 = (N — 
1) _1 J2i( a i ~ ^) 2 ' f° r the variance is used, an indication 
that the unweighted variance in this data is due primarily 
to measurement error, not intrinsic scatter. 



Non-uniform rotation. Now suppose we wish to test 
if there is a single Im spherical-harmonic variation in 
a{0,4>): i.e., that a(9,4>) = ai m Yi m (9, </>), for some 
given I and m. Then each data point would provide 
an estimator, a\ m — (ai — a)/Yi m (9i, (f>i), with variance 
(|cJ; m | 2 ) = af/\Yi m (9i, 4>i)\ 2 - I include the a term in the 
estimator to avoid confusing a higher moment (e.g., a 
dipole) with uniform rotation in case of limited or irregu- 
lar sky coverage. It should become irrelevant in an ideal 
experiment, with N — > oo and a population of sources 
spread uniformly throughout the sky. Note that a\ m is 
complex, and the variances to the real and imaginary 
parts are each (|S?J m | 2 ) /2. 

The minimum-variance estimator a~i m obtained from 
all N data points is obtained by adding all the N in- 
dividual 35 estimators with inverse-variance weighting; 



i.e. 



aim = 



-Y,t 



E 



\Y,r. 



(3) 

with variance given by the inverse of the denominator in 
this expression. The results of such an analysis of the 9 
sources m Table Hare presented in Table HH There is no 
evidence for any nonzero cv; m for I < 2. 
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TABLE II: The measured a irl 
bleU 



obtained from the data in Ta- 



The values of the individual S; m s, for a given I, 
depend on the choice of coordinate system. To test 
for a non-uniform CB in a rotationally-invariant way, 
one must evaluate the rotational invariants Ci = 
S m =-i \airn\ 2 1 [21 + 1). Doing so, no evidence of a 
non- uniform CB is found for the dipole (1 = 1) and 
quadrupole (1 = 2). Rough upper limits to the dipole and 
quadrupole amplitudes can be obtained from the noise: 
y/d/iAn) < 2.3° and ^C 2 /(4ir) < 1.9°. 

Higher-l moments. Since we have in the current anal- 
ysis only 8 data points, it is not really possible to mea- 
sure any ai m s with I > 2. However, if there were 
a nonzero ai m for some high I, it would give rise to 
a scatter in the measured c^s with variance, (a 2 ) = 
(Att)^ 1 J dn [a(9,cj))] 2 = |a 2 m |/(47r). If the sources are 
randomly distributed on the sky at points with angular 
separations A9 3> n/l, then this variance (a 2 ) cannot 
be larger than the variance in the data. The variance 
measured from the data in Table U is roughly (2.9°) 2 , 
but there is a sample error to this variance, of roughly 



3 



y/2/N, where N is the number of (statistically signifi- 
cant) data points, which I estimate to be 5. I therefore 



from which are obtained the rotational invariants, 



take, as a rough upper limit (a 2 ) 1/2 < 3.7°, implying 
c^rnl 2 /4tt < (3.7 ) 2 . The upper limit to any Ci are sim- 
ilar: C//47T < (3.7°) 2 for any individual higher I. 

Stochastic values of a. Theories with a spatially- 
varying a(9, <p) generally predict that a{9, 0) is a real- 
ization of a random field with some given power spec- 
trum C/. Such a theory results in a variance in the 
measured offsets of (a 2 ) = ^2 t (2l + l)Cj/(47r) which, 
again, must be < (3.7°) 2 . For example, suppose some 
theory predicts a scale-invariant spectrum, 1(1 + 1)(7; = 
2Ci exp(— Z 2 // 2 ) from Z = f out to some cutoff moment 
l c with an amplitude parametrized by the dipole mo- 
ment C\. Then approximating for this model (a 2 ) — 
C\ ln(7Z c )/7r, we find the dipole to be constrained to 
C*i/4vr < (1/4) (a 2 ) /ln(7/ c ) ~ (0.7°) 2 / ln(Z c /200). Note 
that the finite angular size 8 of the images limits the ef- 
fective l c < 200 (5/1°) ^ even if the theory allows it to 
be much larger. 

Analysis of individual objects. Early measurements of 
CB |i, 13, 11 1 cons idered simply the offset between the 
image position angle and the mean polarization averaged 
over the entire image. But this averagin g e rases much 
of the information available in the source [12, 13, lal. If 



a high-resolution map of the intensity and polarization 
of a given source is available, then the offset between the 
intensity gradient and the polarization throughout a given 
source can provide a far more precise measurement of the 
mean offset for that particular source. The sensitivity 
to CB from a handful of well-resolved sources can thus 
compete with that of hundreds of unresolved sources. 

Still, one can do better in terms of measuring the offset 
a from a given well-resolved source than prior analyses 
by using techniques developed to quantify temperature- 
polarization correlation functions in the CMB [13, E3], 
and also galaxy-shape correlations induced by weak grav- 
itational lensing [20(. These techniques deal, for exam- 
ple, with the ambiguity in the direction of the linear po- 
larization and also with optimizing low signal-to-noise 
measurements. They allow the full two-point intensity- 
polarization correlations to be used, rather than sim- 
ply the correlations at zero lag (as prior analyses have 
used). They provide additional information on a from 
the polarization alone, even without cross-correlation 
with the intensity, something that could not be done 
with the more heuristic intensity-gradient-polarization 
correlation. These techniques are analogous to those for 
measuring a uniform rotation angle with the CMB y|. 

Suppose that we have a resolved map of the intensity 
P\9 x ,9 y ) and Stokes parameters Q(9 x ,9 y ) and U(9 x ,9 y ) 
of a given radio source, where 9 X and 9 y are coordi- 
nates in the image plane. One first Fourier transforms, 
IQx, ly) = I d 2 9e ltS I{9), and similarly for Q(F) and [7(F), 
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ill 



l 2 y)U(l) 



I 2 



(4) 



From these, the six power spectra P l xx — (X(l)X'(l)*) 
(where the angle brackets denote an ensemble average) 
can be obtained, where {X,X'} = {I,E,B}. The B 
modes have opposite parity from the I and E modes, and 
so we should have P{ B = P BB — 0, if parity is preserved. 
While any given source may in principle have some hand- 
edness, and thus possibly nonzero IB or EB correlations, 
there should be no preference for a given handedness 
when averaging over many sources. More importantly, 
the existing measurements, which show that the offset a 
is small in the sources where it is measured, suggest that 
IB and EB correlations will be small. 

If CB rotates the polarization by an angle a, then 
part of the E mode is rotated into a B mode, 8B{1) = 
E(l)$m2a ~ 2aE(l), thus inducing nonzero Pf B = 
2aP t IE and P EB = 2aP l EE . A rotation-angle estimate 
is then obtained by comparing the measured IB and EB 
correlations with the measured IE and EE correlations, 
respectively. 

A similar analysis can be done, alternatively and equiv- 
alently, using temperature-polarization two-point corre- 
lation functions [13, [13] ■ This involves taking all pairs 
(^1,^2) °f points in the map, and then measuring cor- 
relations between the intensity I and Stokes parame- 
ters Q r and U r measured in a coordinate system that is 
aligned with the line connecting the two points. Again, 
symmetry considerations suggest that, in the absence 
of CB, (J(&)ET r (&)) = (Qr(9i)U r (9 2 )) = 0. If the 
polarizations are rotated by an angle a, then these 
parity-odd correlations are induced, with magnitudes 
(I(0i)EZr(&)) = 2a(Wi)Qr{0 2 )) and (Q r (6i)U r (6 2 )) = 
2a(Q r (Oi)Q r (02)). The decision as to whether to use 
power spectra or correlation functions will depend on the 
noise properties of the map. 

Discussion. Here I have discussed measurements of a 
CB rotation of the linear polarization that varies as a 
function of position on the sky and derived rough upper 
limits to rotation-angle multipole moments and power 
spectra. I discussed how the analysis of future high- 
resolution intensity-polarization images of high-redshift 
sources can be optimized with techniques analogous to 
those in CMB-polarization studies. 

The analysis presented here is meant primarily to be 
illustrative. The existing data are far from optimized for 
this particular measurement. First of all, I used only 9 
sources at redshifts z > 2, and the statistical weight is 
dominated by only half of those. Moreover, they are not 
uniformly spread on the sky (which is why the errors to 



4 



the different m moments for a given I vary so widely), and 
this could give rise to pitfalls. Had my analysis found ev- 
idence for a signal, this may have been cause for concern. 
But given that the results are null, the derived upper 
limits are probably sound. 

Although a comparable sensitivity to a position- 
dependent rotation can in principle be obtained from ex- 
isting CMB data, the analysis is difficult and has not yet 
been done. The simple exercise I have performed here 
is thus the strongest existing constraint to a position- 
dependent rotation, at least for a rotation that occurs at 
redshifts < z < 2. A model that predicts rotation at 
3 < z < 1100 could still produce a signal in the CMB 
without violating the constraint I have derived. Like- 
wise, a slightly stronger constraint can probably be ob- 
tained from the radio-galaxy data in Ref. [IH , although 
for lower-redshift sources, and thus over a smaller base- 
line. Whether that constraint would be competitive with 
the one I have derived would, again, depend on the red- 
shift dependence of a in any given model. Of course, if 
one has a particular model that makes a specific predic- 
tion for the redshift dependence of a(n; z), then the data 
from sources at all redshifts can be combined to provide 
optimal estimators for the parameters of the model. 

Improvements to the illustrative analysis I have done 
here should be straightforward. The error to the ai m 
should scale simply as N^ 1 / 2 with the number N of 
sources, assuming similar image qualities to those ob- 
tained so far. Thus, for example, if the sample I used 
of N ~ 4 well-measured offsets can be improved to 
N ~ 400, the sensitivity will be competitive with that 
expected from Planck. However, progress can be acceler- 
ated, beyond N~ x / 2 , if more precise offset measurements 
can be obtained for at least some of these individual 
sources, either from better images, an improved analysis, 
or both. The goal of identifying the new physics respon- 
sible for cosmic acceleration will hopefully motivate such 
empirical investigations. 
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